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ABSTRACT 

We reexamine nonlinear diffusive shock acceleration (DSA) at cosmological 
shocks in the large scale structure of the Universe, incorporating wave-particle 
interactions that are expected to operate in collisionless shocks. Adopting simple 
phenomenological models for magnetic field amplification (MFA) by cosmic-ray 
(CR) streaming instabilities and Alfvenic drift, we perform kinetic DSA simu- 
lations for a wide range of sonic and Alfvenic Mach numbers and evaluate the 
CR injection fraction and acceleration efficiency. In our DSA model the CR ac- 
celeration efficiency is determined mainly by the sonic Mach number Mg, while 
the MFA factor depends on the Alfvenic Mach number and the degree of shock 
modification by CRs. We show that at strong CR modified shocks, if scattering 
centers drift with an effective Alfven speed in the amplified magnetic field, the 
CR energy spectrum is steepened and the acceleration efficiency is reduced sig- 
nificantly, compared to the cases without such effects. As a result, the postshock 
CR pressure saturates roughly at ~ 20 % of the shock ram pressure for strong 
shocks with Mg > 10. In the test-particle regime {Ms < 3), it is expected that 
the magnetic field is not amplified and the Alfvenic drift effects are insignificant, 
although relevant plasma physical processes at low Mach number shocks remain 
largely uncertain. 



Subject headings: acceleration of particles — cosmic rays — galaxies: clusters: 
general — shock waves 
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INTRODUCTION 



It is expected that hierarchical gravitational clustering of matter induces shock waves in 
baryo nic gas in the large-scale structure (LSS) of the Universe (IKang et al.lll996l : lMiniati et al 
2000|). Simulations for the LSS formation suggest that strong shocks [Mg ^ 10) form in rel- 
atively cooler environments in voids, filaments, and outside cluster virial radii, while weak 
shock {Ms 5, several) are produced by mergers and flow motions in hotter intracluster me- 



dia dCMs) flRvu et al. 



Hoeft et al. 2008: Vazza et al. 



2003 



'frommer et al 



2009 



2006 



Kang et all 120071 ISkillman et al.l 12008 



Briiggen et al.l |2012| ) . Observationally, the existence 



of such weak shocks in ICMs has been revealed through temperature jumps in the X-ray 
emitting gas and Mpc-scale re l ics with radio spect ra softening downstream of the shock (see 
Markevitch fc Vikhliniru 120071 : iFeretti et al.ll2012l . for reviews). These cosmological shocks 
are the primary means through which the gravitational energy released during the LSS for- 
mation is dissipa ted into the gas entropy, magnetic fleld, turbulence and nonthermal particles 
JRvu et al.ll2"008l ). 



In fact, shocks are ubiquitous in astrophysical environments from the heliosphere to 
galaxy clusters and they are thought to be the ma in 'cosmic accelerators' of high energy 



cosmic-ray (CR) particles (IBlandford fc Eichlerlll987l ). In diffusive shock acceleration (DSA) 



theory, suprathermal particles are scattered by magnetohydrodynamic (MHD) waves and 



(Bell 


1978: 


Drurv 


1983: 



are advected downstream, some suprathermal particles energetic enough to swim against 
downstream turbulent waves can cross the shock and be injected into the Fermi flrst-order 
process. Then these streaming CRs generate resonant waves via two-stream instability and 
nonresonant waves via CR current- driven instability, which in turn amphfy turbulent mag- 



netic flelds in the preshock region flBelll Il978l : iLucek fc Belli I2OOOI : iBelll 12004 ISchure et al. 



2OI2I ). Thin X-ray synchrotron emitting rims observed in several young supernova remnants 
(SNRs) indicate that CR electrons are accelerated to 10-100 TeV an d cool radiatively in 



the magnetic field o f several 100 /iG behind the forward shock (e.g., iParizot et al.l 12006 



Reynolds et al.ll2012[ ). This provides clear evidence for efficient magnetic field amplification 
during CR acceleration at strong CR modified shocks. 

These plasma physical processes, i.e., injection of suprathermal particles into the CR 
population, excitation of MHD waves and amplification of turbulent magnetic fields via 
plasma instabilities, and further acceleration of CRs via Fermi first-order process are im- 
portant ingredients of DSA and should op erate at all types of as t rophysical shocks includ- 



Schure et al. 



ing cosmological shocks in the LSS (e .g., iMalkov fc Druryl 1200 ll : IZweibel fc EverettI 12010 



2OI2I : iBriiggen et al.ll2012l ) . In addition, relativistic particles can be accelerated 
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stochastically by MHD turbulen c e, most likely driven in ICM s of merging clusters (jPetrosian 



2001 : Cassano &: Brunetti 2005 : Brunetti &: Lazarian l2007f) . CR s can be also injected mto 



the intergalactic spac e by radio galaxies (Kronberget al.l l200l[ ) and through winds from 
star-forming galaxies (jVolk fc Atoyanlll999l ). and later re-accelerated by turbulence and/or 
shocks. Diffuse synchrotron emission from radio halos and relics in galaxy clusters indi- 
cates the presence of GeV electrons gyrating i n AtG-level magnetic fie l ds on Mpc scales (e.g., 
Carilh fc Tavloi1l2002[ [Covoni &: Ferettilbooi Ivan Weeren et al.lboioiETang et al.ll2012[ ). On 
the other hand, non-detection of 7-ray emission from galaxy clusters by Fermi-LAT and VER- 
ITAS observations, combined with radio halo observations, puts rather strong constraints 
on the CR proton population and the magnetic field strength in ICMs, if one adopts the 
"hadronic" model, in which inelastic collision s of CR protons with ICM protons produce 



the radio emitting electrons and the tt" dec ay f lAckermann et al.l I2OIOI : iDonnert et al.l 12010 



Jeltema fc Profumoll201ll : lArlen et al.ll2012[ ). Alternatively, in the "re-acceleration" model, 
in which those secondary electrons produced by p-p collisions are accelerated further by 
MHD turbulence in ICMs, the CR proton pressure not exceeding a few % of the gas thermal 
pressure could be consistent with both th e Fermi-LAT upper l imits from the GeV 7-ray flux 



and the radio properties of cluster halos (IBrunetti et al. 



20m . 



Recently, amplification of turbulent magnetic fields via plasma instabilities and injec- 
tion of CR protons and electrons at non-relativistic collisonless sh ocks have been studied, 
using Part i cle-in-Cell (PIC ) and hybrid plasma simul ations (e.g. iRiquelme fc Spitkovsky 



20091 . 12OIII : iGuo et al.ll2010l : ICarate fc Spitkovskvl 120121 ). In PIC simulations, the Maxwell's 
equations for electric and magnetic fields are solved along with the equations of motion for 
ions and electrons, so the full wave-particle interactions can be followed from first principles. 
However, extremely wide ranges of length and time scales need to be resolved mainly because 
of the large proton to electron mass ratio. In hybrid simulations, only the ions are treated 
kinetically while the electrons are treated as a neutralizing, massless fluid, alleviating severe 
computational requirements. However, it is still prohibitively expensive to simulate the full 
extent of DSA from the thermal energies of background plasma to the relativistic energies of 
cosmic rays, following the relevant plasma interactions at the same time. So we do not yet 
have full understandings of injection and diffusive scattering of CRs and magnetic field am- 
plification (MFA) to make precise quantitative predictions for DSA. Instead, most of kinetic 
DSA approaches, in which the diffusion-convection equation for the phase-space distribution 
of particles is solved, commonly adopt phenomenological models that may emulate some of 



those processes (e.g. , iKang et al.ll2002l : iBerezhko et al.ll2009t iPtuskin et al.ll2010t iLee et al. 



2OI2I : ICapriolil l2012t iKangi |2012| ). Another approximate method is a steady-state Monte 
Carlo simulation approach, in which parameterized models for particle diffusion, growth of 
self-generated MHD turbulence, wave dissipation and plasma heating are implemented (e.g.. 



-4- 



Vladimirov et al.lboOSf ). 



In our previous studies, we performed DSA simulations of CR protons at cosmological 
shocks, assuming that the magnetic field strength is uniform in space and constant in time, 
and presented the time-asymptotic values of fractional thermalization, S(M,c,), and fractional 
CR acceleration, rjlMg), as a function of the sonic Mach number Mg ( iKang &: JonesI 120071 : 
Kang et al.ll2009l ). These energy dissipation efficiencies were adopted in a post-processing 
step for structur e formation simulatio r is in order to estim ate the C R generation at cosm olog- 
ical shocks (e.g., ISkillman et al.ll2008l : IVazza et al.ll2009l ). Recently, IVazza et al.l ( 12012| ) have 
used those efficiencies to include self-consistently the CR pressure terms in the gasdynamic 
conservation equations for cosmological simulations. In this paper, we revisit the problem 
of DSA efficiency at cosmological shocks, including phenomenological models for MFA and 
drift of scattering centers with Alfven speed in the amplified magnetic field. Amplification 
of turbulent magnetic fields driven by CR streaming in stabilities is included through an ap- 
proximate, analytic model suggested by ICapriolil ( 12012| ). As in our previous works, a thermal 
leakage injection model and a Bohm-like diffusion coefficient {n{p) oc p) are adopted as well. 

This paper is organized as follows. The numerical method and phenomenological models 
for plasma physical processes in DSA theory, and the model parameters for cosmological 
shocks are described in Section 2. We then present the detailed simulation results in Section 
3 and summarize the main conclusion in Section 4. 



DSA MODEL 



In the diffusion approximation, where the pitch-angle distribution of CRs is nearly 
isotropic, the Fokker-Plank equation of the particle distribution function is reduced to the 
following diffusion-convection equation: 

pd{u + u^) df 
3 



dx 



d_ 

dp dx 



■ . .Of 
^^^'^'~dx 



(1) 



where f{x,p,t) is the isotropic part of the pitch- angle averaged CR distribution function, 
K,{x,p) is the spatial diffusion coefficient along the direction parallel to the mean magnetic 
f ield and M,„ is the drift speed of local Alfvenic wave turbulence with respect to the plasma 
( lSkilling]|l975[ ). Here, we consider quasi-parallel shocks in one-dimensional planar geometry, 
in which the mean magnetic field is roughly parallel to the flow direction. The flow velocity, 
u, is calculated by solving the momentum conservation equation with dynamical feedback of 
the CR pressure and self-generated magnetic flelds. 



d{pu) , d{pu'+Pg + P, + PB) 



dt 



dx 



0. 



(2) 
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The CR pressure, Pc, is calculated self-consistently with the CR distribution function /, 
while the magnetic pressure, Pb, is calculated according to our ph enomenological mod el for 
MFA (see Section 2.4) rather than solving the induction equation (jCaprioli et al.ll2009l ). We 
point that the dynamical effects of magnetic field are not important with Pb ^ O.OlpoW^. 
T he details of o ur PS A numerical code, the CRASH (Cosmic- Ray Amr SHock), can be found 



m 



Kang et all fl2002f ). 



2.1. Thermal Leakage Injection 



Injection of protons from the postshock thermal pool into the CR population via wave- 
particle interactions is expected to depend on several properties of the shock, including 
the sonic and Alfvenic Mach numbers, the obliquity angle of mean magnetic field, and the 
strength of pre-existing and self-excited MHD turbulence. As in our previous studies, we 
adopt a simple phenomenological model in which particles above an "effective" injection 
momentum pinj get injected to the CR population: 

Pinj ^ 1.17mpU2 



^B 



(3) 



where = Bq/B± is the ratio of the mean magnetic fiel d along the shock normal, Bn, to the 



ampl itude of the postshock MHD wave turbulence, B± (IMalkov fc Druryll2001t iKang et al. 



2OO2I ). This injection model reflects plasma physical arguments that the particle speed must 
be several times larger than the downstream flow speed, U2, depending on the strength of 
MHD wave turbulence, in order for suprathermal particles to leak upstream across the shock 
transition layer. Since the physical range of the parameter is not tightly constrained, 
we adopt = 0.25 as a canonical value, which results in the injected particle fraction, 
C, = ^1^,2/^2 ~ 10~'^ — 10~^ for Ms > 3 (see Figure 3 below). Previous studies showed 
that DSA saturates for ^ > 10"^, so the acceleration e fficiency obtained here may repre sent 



an upper limit for the efficient injection regime (e.g., iKang et al.l l2002t ICapriolil |2012| ). In 
fact, this injection frac tion is similar to the commonly adopted values for nonlinear DSA 
modeling of SNRs (e.g.. iBerezhko et al.ll2009l ). If we adopt a smaller value of for stronger 



wave turbulence, Pi^j has to be higher, leading to a smaller injection fraction and a lower 
acceleration efficiency. 



2.2. Bohm-like Diffusion Model 



In our model, turbulent MHD waves are self-generated efficiently by plasma instabilities 
driven by CRs streaming upstream in the shock precursor, so we can assume that CR particles 
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are resonantly scattered by Alfven waves with fully saturated spectrum. Then the particle 
diffusion can be approximated by a Bohm-like diff usion coefficient, Kb ~ {l/3)rgV, but with 
flattened non-relativistic momentum dependence (IKang fc JonesI 120071 ): 

* Bo 



mpC^/{3eBo) = (3.13 x 1022cm2s-^)5o-i 



P 



TTlpC 



(4) 



where n* = rripC /(3e5o) = (3.13 x 10^^cm^s-^)5o-\ and Bq is the magnetic field strength far 
upstream expressed in units of /iG. The strength of the parallel component of local magnetic 
field, B\\{x), will be described in the next section. Hereafter, we use the subscripts '0', T, 
and '2' to denote conditions far upstream of the shock, immediate upstream and downstream 
of the subshock, respectively. 



2.3. Magnetic Field Amplification 

It was well known that CRs streaming upstream in the shock precursor excite resonant 
Alfven waves with a wavelength (A) comparable with the CR gyroradius (r^), and turbu- 



lent magnetic fields can be amplified into the nonlinear regime {i.e., 6B ^ Bq) (jBelllll978 
Lucek fc Bedl2000l ). Later, it was discovered that the nonresonant (A ^ r^), fast-growing 
instability driven by the CR current {jcr = ericrUs) can amplify the magnet ic field by orders 



of magnitude, up to the level consistent with the thin X-ray rims at SNRs (|Bel||200J). Sev- 
eral plasma simulations have shown that both B\\/Bq and B \ / Bn can increase by a factor of 
up to ^ 10 — 45 yia the Bell's CR current-driven instability (IRiquelme fc SpitkovskyI 12009 



2010l : lOhira et al.ll2009l ). Moreover, it was suggested that long- wavelength magnetic fluc- 



tuations can grow as well in the p resence of short- sc ale, circula rly-polarized Alfveri waves 



excited by the Bell- type instability (IBykov et al.ll201ll ). Recently, iRogachevskii et al.l (120121 ) 



have also shown that large-scale magnetic fluctuations can grow along the original fleld by 
the a effect driven by the nonresonant instability and both the parallel and perpendicular 
components can be further amplifled. There are several other instabilities that may amplify 



fllamentation, and acoustic instabilities (e.g. 



Beresnvak et al. 


2009: 


Drurv & Downes 


2012 



Schure et al.ll2012l ). Although Bell's (2004) original study assumed parallel background mag- 
netic fleld, it turns out that the non-resonant instability operates for all shocks, regardless 
of the inclination an gle between the shock normal and the mean background magnetic fleld 
JSchure et al.ll2012h . and so the isoptropization of the amplifled magnetic fleld can be a 
reasonable approximation ( Riquelme fc Spitkovsky 2009 : Rogachevskii et aPboid ). 



Here, we adopt the pres cription for MFA due to CR streaming instabilities that was 
suggested by ICapriolil (120121 ). based on the assumption of isotropization of the amplifled 
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magnetic field and the effective Alfven speed in the local, amplified field: 6B'^ /{SnpouD = 
(2/25) (1 - f/^/4)2f/-i-5, where 5B = B - Bq and f/ = (n, - is the flow speed in the 

shock rest frame normalized by the shock speed Ug. In the test-particle regime where the 
flow structure is not modified, the upstream magnetic field is not amplified in this model 
{i.e., U{x) = 1). In the shock precursor (x > Xg, where Xg is the shock position), the MFA 



factor becomes 



6B{x) 



[1 - U{xf'^f 



-°o 



"25^^^^'° f/(,x)3/2 ' 

where Ma q = Us/va,o is the Alfvenic Mach number for the far upstream Alfven speed, and 
va,o = -Bo/-\/4vrpo • This model predicts that MFA in creases with M/i n and the precursor 
strength {i.e., degree of shock modification by CRs) ( IVladimirov et al.ll2008l ). In the case 
of a "moderately modified" shock, in which the immediate preshock speed is Ui ~ 0.8, for 
example, the amplified magnetic pressure increases to 5B\/'&ti ^ 6.6 x 1Q~^pqu1 and the 
amplification factor scales as SBi/Bq 0.12Ma o- We will show in the next section that the 
shock structure is modified only moderately owing to the Alfvenic drift, so the magnetic field 
pressure is less than a few % of the shock ram pressure even at strong shocks {Mg > 10). 

For the highest Mach number model considered here, Mg = 100, the preshock ampli- 
fication factor becomes 6Bi/Bq ^ 100, which is somewhat larger than what was found in 
the plasma simulations for the Bell-type current- driven instability ( iRiquelme &: Spitkovsky 
20091 . 12OIOI ) . Con sidering possible MFA beyond the Bell-type i nstability by other large-scale 



instabilities (e.g. iBykov et al.ll201ll : iRogachevskii et al.ll2012l : ISchure et al.l 120121 ). this level 



of MFA may not be out of reach. Note that this recipe is intended to be a heuristic model 
that may represent qualitatively the MFA process in the shock precursor. 

Assuming that the two perpendicular components of preshock magnetic fields are com- 
pletely isotropized and simply compressed across the subshock, the immediate postshock 
field strength can be estimated by 



B2/B, = ^1/3 + 2/3(p2/pi) 

We note that the MFA model described in equations ([5])-( 
coefficient model given by equation (jlj). 



(6) 



is also used for the diffusion 



2.4. Alfvenic Drift 



Resonant Alfven waves excited by the cosmic ray streaming are pushed by the CR 
press ure gradient (dPr/dx) an d propagate against the underlying flow in the shock precursor 
(e.g. ISkilling] Il975l : iBelll Il978l ). The mean drift speed of scattering centers is commonly 



- 8 - 



approximated as the Alfven speed, i.e., Uw,i{x) ~ +va ~ B{x)/ ^47rp(x), pointing upstream 
away from the shock, where B{x) is the local, amplified magnetic field strength estimated by 
equation ([5]). For isotropic magnetic fields, the parallel component would be roughly B\\ ^ 
B{x) I a/3. But we simply use B[x) for the effective Alfven speed, since the uncertainty in this 
model is probably greater than the factor of -v/S (see Section 3 for a further comment on this 
factor). In the postshock region the Alfvenic t urbulence is probably relatively balanced, so 
the wave drift can be ignored, that is, 2 ~ (|Joneslll993l ). Since the Alfvenic drift reduces 
the velocity difference between upstream and downstream scattering centers, compared to 
that of the bulk flow, the resulting CR spectrum becomes softer than estimated without 
considering the wave drift. Here, we do not consider loss of turbulent magnetic energy and 
gas heating due to wave dissipation in order to avoid introducing additional free parameters 
to the problem. 



2.5. Set-up for DSA Simulations 



Previous studies have shown that the DSA efficiency depends primarily on the shock 
sonic Mach number (IKang et al.l 120071 ). So we considered shocks with a wide range of the 
sonic Mach number, = 1.5 — 100, propagating int o the intergalactic medium (IGM) of 
different temperature phases, Tq = 10^ — 5 x 10^ K (IKang et al.ll2005l ). Then, the shock 
speed is given by = (150 km s-^)M,(To/106K)1/2. 

We specify the background magnetic field strength by setting the so-called plasma beta, 
/3p = Pg/PB, the ratio of the gas pressure to the magnetic pressure. So the upstream mag- 
netic field strength is given as B^ = 8iTPg//3p, where (3p ~ 100 is taken as a canonical value 
in ICMs (see, e.g.. lRyu et al.ll2008l ). Then, the ratio of the background Alfven speed to the 
sound speed, va,o/cs = ^/2J{fip^ (where 7^ is the gas adiabatic index), which determines 
the significance of Alfvenic drift, depends only on the parameter fip. Moreover, the up- 
stream Alfvenic Mach number, Maa = Us/va,o = Ms\f^PlgJ^-, controls the magnetic field 
amplification factor as given in equation For /3p = 100 and 7^ = 5/3, the background 
Afiven speed is about 10 % of the sound speed, i.e., t>^o = 0.1 Ic^ (independent of Mg and 
To), and Ma^q = O.IM^. For a higher value /3p [i.e., weaker magnetic fields), of course, the 
Alfvenic drift effect will be less significant. 

With a fixed value of /3p, the upstream magnetic field strength can be specified by the 
upstream gas pressure, UhoTq, as follow: 



Bq = 0.28 fiG 



103 cm-3K 



1/2 



100 

1^ 



1/2 



(7) 
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We choose the hydrogen number density, nH,o = 10~^ cm~^, as the fiducial value to obtain 
specific values of magnetic field strength shown in Figures 1-2 below. But this choice does 
not affects the time asymptotic results shown in Figures 3-4, since the CR modified shock 
evolves in a self-similar manner and the time-asymptotic states depend primarily on M<, and 
Ma ft, independent of the specific value of Bq. 

Since the tension in the magnetic field lines hinders Bell's CR current- driven insta- 
bility, MFA occurs if the background field strength satisfies the condition, Bq < Bg = 
(0.87 iiG){ncrUsY^'^ (IZweibel &: Everettll201Cll ). For a typical shock speed of Ug ~ 10'^ km s~^ 
formed in the IGM with nH,o ~ 10~^ — 10~^ cm~'^ with the CR injection fraction, ~ 
10~^ — 10"'^, the maximum magnetic field for the growth of nonresonant waves is roughly 
~ 0.1 — 1 yuG. The magnetic field strength estimated by equation ([7]) is i?o ~ 0.28 /iC 
for UHfi = 10"^ cm-3 and Tq = W K (ICMs) and Bq ^ 10'^ fiG for n^.o = 10"^ cm'^ 
and To = 10^ K (voids). Considering the uncertainties in the model and the parameters, 
it seems reasonable to assume t hat MFA via CR strearn ing instabilities can be effective at 
cosmological shocks in the LSS (IZweibel fc Everettll2010l ). 



In the simulations, the diffusion coefficient, n* in equation (jl]), can be normalized with 
a specific value of Kq- Then, the related length and time scales are given as = i^o/us and 
to = t^o/ul, respectively. Since the fiow structure and the CR pressure approach the time- 
asymptotic self-similar states, a specific physical value of Ko matters only in the determination 
of Pmax/T^pC ~ O-lu^t/K* at a given simulation time. For example, with = IO^k*, the 
highest momentum reached at time t becomes Pmax/fnpC ~ 10^ {t/ to). 

It was suggested that non-linear wave damping and dissipation due to ion-neutral colli- 
sions may weaken stochastic sca tterings, leading to slower acc eleration and escape of highest 
energy particles from the shock (jPtuskin fc Zirakashvilill2005l ). Since these processes are not 
well understood in a quantitative way, we do not include wave dissipation in the simulations. 
Instead we implement a free escape boundary (FEB) at an upstream location by setting 
f{xFEB,p) = at xfeb = 0.5 lo, which may emulate the escape of the highest energy parti- 
cles with the diffusion length, k{p)/us > xfeb- Under this FEB condition, the CR spectrum 
and the shock structur e including the precursor approach the time-asymptotic states in the 
time scale of t/to ~ 1 (lKang||2012l ). 



As noted in the introduction, CR protons can be accelerated by merger and accretion 
shocks, injected into the intergalactic space by star forming galaxies and active galaxies, 
and accelerated by turbulence. Because of long life time and slow diffusion, CR protons 
should be accumulated in the LSS over cosmological times. So it seems natural to assume 
that ICMs contains pre-existing populations of CR protons. But their nature is not well 
constrained, except that the pressure of CR protons is less than a few % of the gas thermal 
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pressure (lArlen et al.l 120121 : iBrunetti et al.ll2012[ ). For a model spectrum of pre-existing CR 
protons, we adopt a simple power-law form, fo{p) = fprc-{p/Pm})~^ for p > Pinp with the slope 
s = 4.5, which corresponds to the slope of the test-particle power-law momentum spectrum 
accelerated at M = 3 shocks. We note that the slope of the CR proton spectrum inferred 
from the radio spectral in dex {i.e., a/j ~ (s — 2)/2) of cluster halos ranges 4.5 ^ s < 5 (e.g., 
Jeltema &: Profumdl201ll ). The amplitude, /pre, is set by the ratio of the upstream CR to 



gas pressure, R = Pc,o/Pgfl, where R = 0.05 is chosen as a canonical value. 

Table 1 lists the considered models: the weak shock models with Tq > 10^ K, the 
strong shock models with Tq = 10^ — 10^ K, and the strongest shock models with Tq = 10^ 
K represent shocks formed in hot ICMs, in the warm-hot intergalactic medium (WHIM) 
of filaments, and in voids, respectively. Simulations start with purely gasdynamic shocks 
initially at rest at Xg = 0. 



3. DSA SIMULATION RESULTS 

Figures 1-2 show the spatial profiles of magnetic field strength, B{x), and CR pressure, 
Pc{x), and the distribution function of CRs at the shock location, gs{p), at t/to = 0.5, 1,2 
for models without or with pre-existing CRs: from top to bottom panels, Mg = 3 and 
To = 5 X lO^K, M, = 5 and Tq = lO^K, M, = 10 and Tq = lO^K, and M, = 100 and 
To = lO^K. Note that the models with Mg = 3 — 5 represent shocks formed in hot ICMs, 
while those with Mg = 10 and 100 reside in filaments and voids, respectively. 

The background magnetic field strength corresponds to Bq = 0.63, 0.28, 0.089, and 8.9 x 
10"'^ fiG for the models with Mg = 3, 5, 10, and 100, respectively, for the fiducial value of 
^H.o = 10~^ cm~'^ (see equation ([7])). With our MFA model the postshock field can increase 
to i?2 ~ 2 — 3 /xG for all these models, which is similar to the field strengths observed in radio 
halos and radio relics. The postshock CR pressure increases with the sonic Mach number, 
but saturates at Pc,2/{pou1) ~ 0.2 for Mg > 10. One can see that the precursor profile and 
gs{p) have reached the time-asymptotic states for t/to ^ 1 for the Mg = 100 model, while 
the lower Mach number models are still approaching to steady state at t/to = 2. This is 
because in the Mg = 100 model, by the time t/to ^ 1 the CR spectrum has extended to Pmax 
that satisfies the FEB condition. For strong shocks of Mg = 10 — 100, the power-law index, 
q = —dlnf/dlnp, is about 4.3 - 4.4 at p ~ m^c instead of g = 4, because the Alfvenic drift 
steepens the CR spectrum. 

For the models with pre-existing CRs in Figure 2, the pre-existing population is im- 
portant only for weak shocks with Mg < 5, because the injected population dominates in 
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shocks with higher sonic Mach numbers. As mentioned in the Introduction, the signatures 
of shocks observed in I CMs through X-ray and radio observations can b e interpreted by low 
Mach number shocks (IMarkevitch fc Vikhhniru 120071 : iFeretti et aLll2012l ). In particular, the 
presen ce of pre-existing CRs is expect to be crucial in explaining the observations of radio 
rehcs kang et al.1l2012l ). 



Figure 3 shows time-asymptotic values of downstream gas pressure, Pg^2, and CR pres- 
sure, Pc,25 in units of pou^, density compression ratios, ai = pi/po and (T2 = P2/ Po, the 
ratios of amplified magnetic field strengths to background strength, -B2/-B0 and Bi/Bq, and 
postshock CR number fraction, = Tier, 2/^2, as a function of Mg for all the models listed in 
Table 1. We note that for the models without pre-existing CRs (left column) two different 
values of To (and so Ug) are considered for each of Mg = 3, 4, 5, 10, 30, and 50 models, in 
order to explore the dependence on To for a given sonic Mach number. The figure demon- 
strates that the DSA efficiency and the MFA factor are determined primarily by Mg and 
Maa, respectively, almost independent of Tq. For instance, the two Mach 10 models with 
To = lO^K (open triangle) and lO^K (filled triangle) show the similar results as shown in the 
left column of Figure 3. But note that the curves for Pcr-,2 and ^ increase somewhat unevenly 
near Mg 4 — 7 for the models with pre-existing CRs in the right column, because of the 
change in To (see Table 1). 

At weak shocks with Mg < 3, the injection fraction is < 10~^ and the CR pressure 
is Pc,2/ PquI ^ 5 X 10~^ without pre-existing CRs, while both values depend on Pc,i in the 
presence of pre-existing CRs. Since the magnetic field is not ar nplified in the test-p article 
regime, these results remain similar to what we reported earlier in lKang fc Ryul (120111 ). For a 
larger va lue of e^, the injection fraction and the CR acceleration efficiency would increase. As 
shown in lKang &: JonesI (l2007l ). however, ^ and Pcr,2/ Po^l depend sensitively on the injection 
parameter for Mg < 5, while such dependence becomes weak for Mg > 10. Furthermore, 
there are large uncertainties in the thermal leakage injection model especially at weak shocks. 
Thus it is not possible nor meaningful to discuss the quantitative dependence of these results 
on e^, until we obtain more realistic pictures of the wave-particle interactions through PIC 
or hybrid plasma simulations of weak collisionless shocks. 

In the limit of large M^, the postshock CR pressure saturates at Pc,2 ~ 0.2po^s' the 
postshock density compression ratio at (T2 ~ 5, and the postshock CR number fraction at 
^ ^ 2 X 10-3. The MFA factors ar e Pi/Pn ~ O . I2M4 n ~ Mg and P2/P0 ~ 3M, for Mg > 5, 
as expected from equation ([5]). In iKang et al.l (120071 ) we found that Pc,2 ~ 0.55poW^ in the 



limit of large Mg, when the magnetic field strength was assumed to be uniform in space and 
constant in time. Here we argue that MFA and Alfvenic drift in the amplified magnetic field 
steepen the CR spectrum and reduce the DSA efficiency drastically. 



-VI- 



Again, the presence of pre-existing CRs (right column) enhances the injection fraction 
and acceleration efficiency at weak shocks of Mg < 5, while it does not affect the results at 
stronger shocks. Since the upstream CR pressure is Pc,o = O.OSPg^o = (0.03/Ms)po^is these 
models, the enhancement factor, Pc,2/-fc,o ~ 1-5 — 6 for Mg < 3. So the DSA acceleration 
efficiency exceeds only slightly the adiabatic compression factor, a'^", where 7c ~ 4/3 is the 
adiabatic index of the CR population. 



As in iKang et al.l (120071 ). the gas thermalization and CR acceleration efficiencies are 
defined as the ratios of the gas thermal and CR energy fluxes to the shock kinetic energy 
flux: 

^ ' (l/2)po«f ' (l/2)po«f ^ ' 

where Cg and are the gas thermal and CR energy densities. The second terms inside the 
brackets subtract the effect of adiabatic compression occurred at the shock. Alternatively, 
the energy dissipation efficiencies not excluding the effect of adiabatic compression across 
the shock can be defined as: 

which may provide more direct measures of the energy generation at the shock. Note that 
7] = rj' for the models with Pc,o = 0. 

Figure 4 shows these dissipation efficiencies for all the models listed in Table 1. Again, 
the CR acceleration efficiency saturates at ?7 ~ 0.2 for M,^ > 10, which i s much lower tha n 



what we reported in the previous studies without MFA (IRyu et al.l 120031 : iKang et al.l 120071 ). 
The CR acceleration efficiency is < 0.01 for weak shocks {Mg < 3) if there is no pre-existing 
CRs. But the efficiency t]' can be as high as 0.1 even for these weak shocks, depending on 
the amount of pre-existing CRs. The efficiency rj for weak shocks is not affected by the new 
models of MFA and Alfvenic drift, since the magnetic field is not amplified in the test-particle 
regime. 

If we choose a smaller value of (3p, the ratio f^.o/cs is larger, leading to less efficient 
acceleration due to the stronger Afvenic drift effects. For example, for /3p ~ 1 [i.e., equipar- 
tition fields), which is relevant for the interstellar mediur n in galaxie s, the CR acceleration 



efficiency in the strong shock limit reduces to ?7 ~ 0.12 (jKangJ |2012|) . On the other hand, 
if we were to choose a smaller wave drift speed, the CR efficiency rj will increase slightly. 
For example, if we choose u^, ~ O.Sva instead of Uy^ ~ va, the value of rj in the high Mach 
number limit would increase to ~ 0.25 for the models considered here. 

On the other hand, if we choose a smaller injection parameter, for example, = 0.23, 
the injection fraction reduces from ,^ = 2.1 x 10""^ to 6.2 x 10~^ and the postshock CR pressure 
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decreases from Pc,2/po^^s = 0.076 to 0.043 for the = 5 model, while ^ = 2.2 x 10"^ to 
3.3 X 10"^ and Pc,2/ Pqu^ = 0.18 to 0.14 for the = 50 model. Considering that the CR 
injection fraction obtained in these simulations > 10~^) is in the saturation limit of DSA, 
the CR acceleration efficiency, t], for M > 10 in Figure 4 should be regarded as an upper 
limit. 



4. SUMMARY 

We revisited the nonlinear DSA of CR protons at cosmological shocks in the LSS, 
incorporating some phenomenological models for MFA due to CR streaming instabilities 
and Alfvenic drift in the shock precursor. Our DSA simulation code, CRASH, adopts the 
Bohm-like diffusion and thermal leakage injection of suprathermal particles into the CR 
population. 

A wide range of preshock temperature, 10^ ^ Tq < 5 x lO^K, is considered to represent 
shocks that form in clusters of galaxies, filaments, and voids. We found that the DSA 
efficiency is determined mainly by the sonic Mach number Mg, but almost independent of 
Tq. We assumed the background intergalactic magnetic field strength, Bq, that corresponds 
to the plasma beta /3p = 100. This is translated to the ratio of the Alfven speed in the 
background magnetic field to the preshock sound speed, va,o/cs = \/6/5l3p ~ 0.11. Then 
the Alfvenic Mach number Mj^ q = ^5/3p/6 Mg determines the extent of MFA (i.e., Bi/Bq), 
which in turn controls the significance of Alfvenic drift in DSA. Although the preshock 
density is set to be nH,o = 10~^ cm~'^ just to give a characteristic scale to the magnetic 
field strength in the IGM, our results for the CR proton acceleration, such as the dissipation 
efficiencies, do not depend on a specific choice of uha- If one is interested in CR electrons, 
which are affected by synchrotron and inverse Compton cooling, the electron energy spectrum 
should depend on the field strength Bq and so on the value of n^.o^b (see equation ([7])). 

The main results of this study can be summarized as follows: 

1) With our phenomenological models for DSA, the injected fraction of CR particles is 
^ ^ 10~^ — 10~^ and the postshock CR pressure becomes 10"^ < -Pc,2/(po'"s) ^ 0.2 for 3 < 
Mg < 100, if there are no pre-existing CRs. A population of pre-existing CRs provides seed 
particles to the Fermi process, so the injection fraction and acceleration efficiency increase 
with the amount of pre-existing CRs at weak shocks. But the presence of pre-existing CRs 
does not affect ^ nor Pc,2 for strong shocks with Mg > 10, in which the freshly injected 
particles dominate over the re-accelerated ones. 

2) The nonlinear stage of MFA via plasma instabilities at collisioness shocks is not fully 
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un derstood yet. S o we adopted a model for MFA via CR streaming instabilities suggested 
by ICapriolil (120 121 ) . We argue that the CR current, jcr ~ G(,o'2nH,oUs, is high enough to 



overcome the magnetic field tension, so the Bell-type instability can amplify tu rbulent mag- 
netic fields at cosmological shocks considered here ( Zweibel fc EverettI |2010[ ). For shocks 



with M > 5, DSA is efficient enough to develop a significant shock precursor due to the 
CR feedback, and the amplified magnetic field strength in the upstream region scales as 
Bi/Bo ^ 0.12Ma,o ^ {(^p/lOOy/'^Ms. This MFA model predicts that the postshock mag- 
netic field strength becomes i?2 ~ 2 — 3 /xG for the shock models considered here (see Table 

!)• 

3) This study demonstrates that if scattering centers drift with the effective Alfven 
speed in the local, amplified magnetic field, the CR energy spectrum can be steepened and 
the acceleration efficiency is reduced significantly, compared to the cases without MFA. As 
a result, the CR acceleration efficiency saturates at 77 = 2ec,r/poul ~ 0.2 for Ms ^10, which 
is sig nificantly lower than what we reported in our previous study, 77 ~ 0.55 (IKang et al. 
20071 ). We note that the value rj at the strong shock limit can vary by ~ 10 %, depending 



on the model parameters such as the injection parameter, plasma beta and wave drift speed. 
Inclusion of wave dissipation (not considered here) will also affect the extent of MFA and 
the acceleration efficiency. This tells us that detailed understandings of plasma physical 
processes are crucial to the study of DSA at astrophysical collisionless shocks. 

4) At weak shocks in the test-particle regime {Mg < 3), the CR pressure is not dy- 
namically important enough to generate significant MHD waves, so the magnetic field is not 
amplified and the Alfvenic drift effects are irrelevant. 

5) Finally, we note that the CR injection and the C R streaming instabilities ar e found 
to be less efficient at quasi-perpendicular shocks (e.g. iGarate &: Spitkovskvl I2OI2I 1 It is 



recognized, however, streaming of CRs is facilitated through locally parallel inclination of 
turbulent magnetic fields at the shock surface, so the CR injection can be effective even 



(Giacalone 


2005; 


Zank et al. 


2006) 



diffus ion coefficient is smaller due to drift motion of particles along the shock surface f lJokipii 



19871 ). In fact, the diffusion convection equation ([T]) should be valid for quasi-perpendicular 
shocks as long as there exists strong MHD turbulence sufficient enough to keep the pitch 
angle distribution of particles isotropic. In that case, the time-asymptotic states of the CR 
shocks should remain the same even for much smaller k{x,p), as mentioned in Section 2.5. 
In addition, the perpe ndicular current-driven instab i hty is found to be effective at quasi- 



perpendicular shocks f lRiquelme fc SpitkovskyI I2OIOI : ISchure et al.ll2012[ ). Thus we expect 



that the overall conclusions drawn from this study should be applicable to all non-relativistic 
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shocks, regardless of the magnetic field inclination angle, although our quantitative estimates 
for the CR injection and acceleration efficiencies may not be generalized to oblique shocks 
with certainty. 
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Table 1. Models 



To(K) P,,o = Pe,o = 0.05P,,o 



5 X 10^ 


M, = 


1.5, 2, 3, 4 


Ms = 


1.5, 2, 3 


10^ 


M, = 


3, 4,5 


Ms = 


4,5 


10^ 


= 


5, 7, 10 


Ms = 


7, 10 


10^ 


Ms = 


10, 20, 30, 50 


Ms = 


20, 30, 50 


10^ 


Ms = 


20, 50, 100 


Ms = 


100 
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, b(mg) , iog(pypnu3 
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Fig. 1. — Magnetic field strength, CR pressure profile and the CR distribution at the shock 
location att/to = 0.5 (dotted lines), 1 (dashed), and 2 (solid) for models without pre-existing 
CRs. See Table 1 for the model parameters. 
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-0.5 0.5 -0.5 0.5 -2 2 4 6 



x/x„ x/x„ log(p/mpC) 

Fig. 2. — Same as Figure 1 except that models with pre-existing CRs are shown. The pre- 
existing population has a power-law spectrum, fp oc that corresponds to an upstream 
CR pressure, Pc,o = O.OSP^^o- In ths right column, the long-dashed lines show the pre-existing 
population. 
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P,o=0 P,,,=0.05Pg,„, s=4.5 




Fig. 3. — Time-asymptotic values of downstream gas pressure, Pg,2, and CR pressure, Pc,2, in 
units of PqU^, density compression ratios, o"i = pi/po and (T2 = P2/P0) the ratios of amplified 
magnetic field strengths to background strength, B2/B0 and Bi/Bq, and postshock CR 
number fraction, ^, as a function of the sonic Mach number, Mg. The left column shows 
the cases without pre-existing CRs, while the right column shows the cases with pre-existing 
CRs. Filled circles are used for the models with Tq = 5 x lO^K, open circles for Tq = lO^K, 
filled triangles for Tq = lO^K, open triangles for Tq = lO^K, and stars for Tq = lO^K. 
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Fig. 4. — Shock dissipation efficiencies in the form of gas and CR energies, 5 and rj (sohd 
hnes), respectively, in equation and 6' and rj' (dotted hnes) in equation OH]) as a function 
of the shock sonic Mach number. For the models with different preshock temperature, Tq, 
the same type of symbols are used as in Figure 3. 



